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ABSTRACT 




Both integrator and canceller types of sampled 
analog comb filters implemented by a canonical second 
order recursive circuit have been studied 
experimentally and theoretically. Two 96 stage CTD 
(charge transfer device) delay lines are used in the 
experimental study. Theoretically, an analysis based 
on a modification of the digital recursive filter 
theory was developed and accounted for the major 
features of the experimental results. 

Three comb filter applications have been 
demonstrated. The first two used the integrator type 
comb filter for the improvement of the signal-co-noise 
ratio of a periodic signal contaminated by noise and 
for the sorting of a periodic signal of selected PRF 
from other periodic signals. The third application 
used a canceller type comb filter to eliminate a 
sinusoidal interference from a desirable signal. 
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DEFINITION OF TERMS 



Several equations are presented throughout this 
thesis. Many of the terms used in these equations are 
generally represented in most texts and papers as Greek 
symbols. To simplify the typing requirements, these terms 
have been redefined. The following list shows the standard 
symbols which are used throughout this text. 



PI = 3.141592654 
f = frequency = 1/T 
w = 2*PI*f 
e = exp 
s = j w 
Z = exp(sT) 

6 

MHz = Mega Hertz = 10 Hertz 

3 

K = Kilo Hertz = 10 Hertz 
db = decibels 

-12 

Pf = Pico Farrads = 10 Farrads 
t = Tao 

t = Any Arbitrary starting Time 
0 

CTD = Charge Transfer Device 

SNR = Signal To Noise ratio 

CZT = chirp-Z-transform 

Transfer Function = V /V 

output input 
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I 



I. 



IN TROD UCTION 



A. GENERAL 



One of the major considerations in system design is 
signal processing. For all types of systems, it is 

necessary to consider how the signal will be represented 
(analog, digital, or sampled analog) , received and 
conditioned to minimize errors and noise; and finally, how 
the desired signal will be recognized and/or discriminated. 
All these factors must be considered whether the system 
under concern is a communication, radar, sonar, accoustical, 
optical, or some other type. Although the signal processing 
techniques vary greatly in all the previously mentioned 
systems, the operations performed have a great deal in 
common. The many processing operations can be grouped into 
four categories; 

Two types of filters; 

1. Nonrecursive 

a. convolvers 

b. Correlators 

2. Recursive 

Two types of spectral analysis; 

1 . Fourier type of transforms 

2. Other types of transforms 
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This thesis is concerned with the sampled analog 
recursive filter. Specifically, the recursive comb filter 
is studied both experimentally and theoretically. The 
feasibility of this filter's application is also 
demonstrated . 

B. RECURSIVE COMB FILTERS 



A comb filter is a class of filters whose transfer 
function is periodic in the frequency domain. There are two 
classifications of comb filters: canceller or bandstop, and 
the integrator or bandpass. The canceller is characterized 
as having a very narrow but a high attenuation region 
between neighboring teeth and good transmission within the 
passband region. The integrator is the opposite of this in 
that the passband region exhibits good signal passage in 
periodically spaced narrow bands of frequencies and a 
relatively wide attenuation band in between neighboring 
teeth. Comb filters have been developed for both analog and 
digital filtering systems. The development of these 
circuits has been based on the feedforward, frequency 
sampling, fast Fourier transform, and recursion technique of 
implementation. The common feature of circuit 
implementaticn for all these devices is that they all 
require delay lines as the basis of their operation. 

The delay devices which are used vary with the systems 
encountered. For example, shift registers are used in 
digital systems whereas sampled analog systems use the 
bucket brigade device (B3D) , charge coupled device (CCD) , or 
the serial analog delay (SAD) device. Due to the sampling 
technique of these sampled analog devices, they exhibit both 
sampled analog properties and some of the digital properties 
as well. 
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The comb filter was analyzed for a sampled analog filter 
to determine its frequency response characteristics. The 
technique used in designing this type of filter was the 
Bilinear Z transform which is discussed in Chapter 2. This 
filter was iuplemented as a recursive filter which is a 
filter that is characterized as a device whose output 
depends on the present input, past output, and its present 
output. It can also be implemented as a non- recursive 
feedforward filter and used for spectral analysis through 
the chir p“Z“transf orm (CZT) . 



C. PURPOSE OF THE STUDY 



During the past two years, recursive comb filters have 

been studied at the Haval Postgraduate school. V. 

[ 10 ] 

lamsa-ad first studied the use of charge transfer device 

delay lines for the implementation of sampled analog 

[7] 

recursive filters. B. Fruend studied the first order 

recursive filter and emphasized the comb feature. L. 

[8 ] 

Saetre extended the study to the second order recursive 

comb filter and demonstrated the application of the 

[6 ] 

canceller type comb filter to MTI radar. S. Holmes 

developed a computer aided design procedure for the 

recursive comb filters and studied their sensitivity to 
device errors and variations. The status of the recursive 
comb filter study before the beginning of this thesis can be 
summarized as follows: 

-Theory- 

The digital recursive filter theory was modified to analyze 
the sampled analog recursive comb filters. While the 
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modified theory accounted for the comb feature of the 
filter, it was not able to explain three major features of 
sampled analog recursive comb filters which are different 
from an ideal comb filter. 

♦The peak of the frequency characteristics do not remain 
constant as in an ideal comb filter. They become frequency 
dependent in the sampled analog implementation. 

♦The frequency characteristics do not remain zero between 
neighboring teeth as in the ideal case. Deterioration of 
these nulls generally increases with frequency. 

♦The shape of the comb teeth do not remain constant as in 
the ideal comb filter. These also change with frequency. 

-experimental study and Application- 

Only the canceller type comb filter was experimentally 
studied in detail and its application to MTI radar was 
demonstrated. The Integrator comb filter was not studied 
thoroughly. 

Therefore, there are three objectives for this thesis; 

1. Improve the theoretical analysis. 

2. Experimentally study the integrator comb filter and 
demonstrate its applications. 

3 . Demonstrate other applications of the canceller comb 
filter and discuss some of its applications. 



D. OOTLINE 



Chapter 1 presents the background and purpose of this 
thesis. Chapter 2 summarizes the basic ideas about digital 
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filters together with a discussion of the mapping transforms 
(H (s) into H (z) ) which ire used for the design of recursive 
comb filters. In Chapter 3, the design procedures for both 
the canceller type and th$ integrator type recursive comb 
filters are discussed. The experimental results with comb 
filters are presented and compared with the theoretical 
analysis in Chapter 4. Chapters 5 , 6 , and 7 present 
applications of comb filters. Chapter 5 specifically 
demonstrates the use of the integrator comb filter as a 
means to improve signal- to-noise ratio of a pulse train 
signal which has been contaminated with white Gaussian 
noise. The integrator comb filter is used as a pulse train 
sorter to select a particular pulse train from among varied 
frequency pulse trains, A third application of the comb 
filter is presented for a canceller comb filter in Chapter 7 
where the device is used to eliminate sinusoidal 
interference from a second, desired sinusoid. Finally, 
Chapter 8 summarizes the results and offers suggestions for 
further studies. 
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ri. THE ORY 



A. INTRODUCTION 



The theory of digital recursive filters has been well 
developed in several text boohs and published papers as well 
as some past theses written hare at the Naval Postgraduate 
School(i.e. see referencss 6, 7, and 8). For this reason the 
theory of recursive filters will not be covered in detail, 
although a general overview of sampled analog theory will be 
presented beginning with the transform theory. This will 
lead into a discussion of the transforms which are used for 
designing both digital and sampled analog filters. 



/ 

B. BASIC TRANSFORM THEORY 



The analysis of analog circuits is highly dependent on 
the Laplace Transform to ease the computations which are 
involved as well as to simplify the circuit synthesis of the 
desired equations. For these same reasons transforms have 
been developed for use in digital circuits. It has been 
found that these same transforms can be used in sampled 
analog circuits. The Z Transforms which have been developed 
can be considered as steppingstones from either the Laplace 
3-plane domain, the discrete time domain (t=nT) , or the 
continuous time domain into the Z-plane domain (see figure 
2 . 1 ) . 
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1 




Figure 2.1 Transform Steppingstones 



C. TRANSFORMS 



There are a number of transforms which can be used to 
change a transfer function in the S-domain to a transfer 
function in the Z domain: 

1. Impulse Invariant Transformation (Standard Z) 

2. Bilinear Z Transformation 

3. Mapping of Differentials 

4. Matched Z Transform 

5 . Others 

The most popular of these are the Impulse Invariant and the 
Bilinear Z transforms shown in figure 2.1. These mapping 
functions are given below: 
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The Standard Z Transform: 



Z=exp (sT) 



The Bilinear Z Transform: 



2( 1 “ Z-^ ) 
T{ 1 + Z“i ) 



D. TRANSFOSM PROPERTIES 



The inherent properties of these types of transforms are 
similar although they do have some pronounced differences 
which make one technique more desirable than the other for 
some applications. These properties of both transforms are 
shown in Tables 1 and 2. 

This thesis utilizes the Bilinear Z Transform technique 
because it maps the entire left half S-plane into the entire 
unit circle and there is no aliasing effect. These were the 
two major reasons for using this transform, other reasons 
are shown in Tables 1 and 2. 
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TABLE 1 



TRANSFORM 


PROPERTIES 


Bilinear Z 


1. Simple algebraic mapping between 

continuous and digital filters. 

2. Maps entire jw axis in the S-plane to 

the unit circle in the Z plane. 

3. Realizable, stable continuous circuits 

map into realizable, stable digital 
systems. 

4. Wideband sharp cutoff continuous filter 

can be mapped into wideband sharp 
cutoff digital filter without 
aliasing. 

5. There is a non-linear relation between 

the analog and digital frequencies 
resulting in frequency warping. 

6. The amplitude response of the analog 

systems being transformed must be 
piece-wise constant, otherwise the 
digital frequency response will be 
warped. (See ref. 3, pg. 222-223) 

7. Neither the impulse nor the phase 

response of the analog filter is 
retained in a digital filter. 
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TABLE 2 



TRANSFORM 


PROPERTIES 


Standard Z 


• 

1. Encounters an inherent aliasing problem 

unless the sampling frequency is high. 

2. The input signal must be band limited 

between -PI/T Hz and PI/T Hz for the 
frequency responses of the analog 
filter and its digital filter design 
equivalent to correspond. 

3. The impulse response of the digital 

filter is a sampled version of its 
analog counterpart. 

4. To use this transformation, the order 

of the numerator must be less than the 
order of the denominator of the analog 
filter to avoid intolerable aliasing 
problems in the digital filter. 
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E. 



PHYSICAL SIGNIFICANCE OF Z 



-1 



-1 

The operation of Z is an important property in 

-1 

designing a sampled data circuit. Mathematically, Z is 

simply a complex variable; but, physically, it represents a 
unit of delay for a sampled data circuit. How then can such 
a function be implemented? There are two techniques which 
are used to implement this delay physically. A digital 
circuit uses a shift register or memory. For sampled analog 
circuits the family of Charge Transfer Devices (CTD) delay 
lines are used. 



F. DIGITAL FILTERS VERSUS SAMPL ED- ANALOG FILTERS 



There are major differences between digital filters and 
sampled-analog filters. However, much of the digital filter 
design theory is carried over into sampled analog filters. 
A digital filter receives an analog signal, quantizes it 
(i.e. represents the signal amplitude by a series of pulses 
which make up a coded binary number corresponding to the 
analog amplitude) and then performs its filtering action. 
This quantizing action is accomplished by an 
analog-to-digital (A-D) converter. It is generally an 
operation which introduces errors and extraneous electrical 
noise into the filter system. A sampled-analog filter 
eliminates this A-D conversion by sampling the signal's 
amplitude and storing a quantity of charge proportional to 
the signal amplitude. The storage (or delay action) is 
accomplished by the charge transfer devices for analog 
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filters and by shift registers for digital filters. 

The relationship between digital filters and sampled- 
analog filters that are interchangeable is seen in the 
design techniques. The well developed digital filter design 
theory is applicable to sampled- analog design as well with 
some modification. 

G. THE UNIT CIRCLE 



When mapping from the S-plane to the Z-plane, it is 
important to consider which regions in the S-plane and 
2-plane will correspond. This knowledge is important to 
consider the stability of the filter designed according to a 
transfer function in the 2 domain. 

For the Bilinear 2 transform, the entire left half 
S-plane maps into the unit circle in the 2-plane, shown in 
figure 2.2, whereas the right half S-plane maps outside the 
unit circle. 




FIGURE 2.2 Bilinear 2 Mapping 
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This is an important result, since the mapping of the 
entire jw axis onto the perimeter of a unit circle preserves 
the frequency selective properties of the analog system. 
The mapping of the left half S-plane inside the unit circle 
in the Z-plane ensures that stable analog filters will map 
into stable discrete time filters. 

The desirable mapping characteristics described above 
are not inherent with the Standard Z transform, shown in 
figure 2.3. The Standard Z transform only maps periodic 
strips of the left half S-plane into the unit circle in the 
Z plane. The right half of the same strip maps outside of 
the unit circle. This is the reason for the possibility of 
aliasing at improper sampling frequencies. These strips are 
2PI/T Hz in width and the first desirable mapping strip is 
between -jPI/T and jPI/T radians in the S-plane. 



/N Imaginary 




Z-Plane 



Figure 2.3 Standard Z Mapping 
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III. FIL T ER DKIGN 



k. SAMPLED ANALOG FILTER DESIGN 



^ • Introducti on 



The design and implementation of a filter can be 
approached from many aspects. Basically, the numerous 
techniques which are used fall into one of two catagories; 

the Direct Method or the Indirect Method. Rabiner and 
[3 ] 

Gold in their excellent textbook ha^^e discussed these 

methods in detail. The procedures used in this thesis were 
a version of the indirect method using the previously 
discussed transforms. Figure 3.1 shows the procedural flow 
of these design techniques. In this section, the design and 
characteristics of the recursive comb filter will be 
covered. 
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DIRECT TECHNIQUE 




a 

MODEL THE ABOVE 
ANALOG FREQUENCY 
DOMAIN FILTER CHARACTERISTIC 
DIRECTLY INTO THE 
Z-DOMAIN 

fon 



INDIREQ TECHNIQUE 



SPECIFICATION IN THE 
FREQUENCY DOMAIN 









APPROXIMATION IN THE 
S-DOMAIN H(s) 
USING THE FOLLOWING 
FILTER PROTOTYPES: 
*BUnERWORTH 
*CHEBYSHEV 
•ELLIPTIC 
•BESSEL 



SAMPLED ANALOG 
FILTER 



a 

VERY SIMILAR 
TO DIGITAL FILTERS 



Implementing 



DIGITAL FILTER ACTIVE FILTERS 

H(s) H(z) VARIOUS CIRCUIT 

STANDARD-Z CONFIGURATION USING 
BILINEAR-Z OP. AMPS. 

MATCHED-Z 

ETC, 

Canonical Realization: 



Hardware : 
Ret I CON 
BBD 
CCL 



Parallel^ Cascade 




Figure 3.1 Design Methods 
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2 . 



Recursive Filters 



A filter which depends on its present input, past 
output, and also its present output is defined as a 
Recursive Filter. This definition implies that there are at 
least two operations which must be performed to synthesize 
such a device; feedback and delay. In addition to these two 
operations, a Recursive Filter also includes multiplication 
and addition. All these operations are evident by studying 
the Z transfer function as seen in equation 1 below: 



a^ + a, Z“i a .2 Z ~2 
1 + bjZ-i + bjZ-z 



1 



This equation is developed from either the analog 
low pass or the analog high pass transfer function using the 
Bilinear Z transformation technique. The derivation of this 
equation is a common one and is found in most text books. 
However, for completeness the reader is encouraged to read 
Appendix B of reference 7 pages 99 through 120. Freund, B. 
[7] 

R. in his thesis has derived, using the Bilinear Z 

transform method, the transfer function for the first and 
second order low and high pass filter. Since the derivation 
of equation 1 has been thoroughly explained by Fruend the 
design techniques here were used to determine the 
coefficients for the proper response that was desired. 

The relationship between the numerator coefficients 
is an important one since the numerator polynomial 
determines the Location of the zeros of the transfer 
response. In the derivation of H(z) the normalized transfer 
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9 



function sho"wed that the relationship between a , a 

0 1 



and a 

2 



is 2 , respectively. As shown below, this gives a 
double zero on the unit circle either at Z=1+j0 or Z=1-j0 
depending upon the sign of the a coefficient. 



For a =+2; 
1 



-1 -2 

U2Z +Z = 0 



2 

Multiply by Z 



2 

Z +2Z+1 = 0 



or 



(Z +1) (Z +1) 
1 2 



0 



therefore 



Z 

1 



= - 1 , 



Z 

2 



= -1 



Similarily for a 



2 



2 

Z -2Z+1 



0 



(Z-1) (Z-1) = 0 

Z = 1, Z =1 
1 2 

Plotting these values of Z on the unit circle reveals that 
there is a double zero at Z= +1 or at Z = -1. This 
indicates that either a low or high pass filter is 
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determined simply by the sign of (plus for low 

minus for high pass). See figure 3.2. Other 

relationships between a , a , and a will result 

0 1 2 

zeros becoming complex, and/or shifting to either 

outside the unit circle resulting in deterioration 
transfer function stop bands* 



pass and 
numerical 
in the 
inside or 
of the 




a. Low Pass Filter 
3 . I = + ^. 0 



FIGORE 3.2 Zero location 




b. High Pass Filter 

a^ = -i.O 

- -f i. 

for low or high pass filter. 



The ”b” coefficients also effect the transfer 
characteristics, although much more differently than the a” 
coefficients. The denominator coefficients affect the shape 
of the pass band region width. By either increasing or 
decreasing the pass band, the 3db cutoff frequency is also 
changed. This change in the pass band shape is accomplished 
by the shifting of the pole locations within the unit circle 
as seen by the calculation below: 
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Denominator equation 



-1 -2 

Ub Z +b Z = 0 
1 2 



Multiplying through by Z 



2 

Z +b Z+b =0 

1 2 



Factoring 



(Z 



1 





2 



Thus the problem of the filter design becomes a 
problem of selecting the proper value of the denominator 
coefficients in the Z transfer function. 



3 . Estimation of th e 3 C oeffici ents 



Insight in determining what the coefficients should 
be for a particular response can be obtained from the 
comparison of the unit circle pole-zero plot and the desired 
- frequency response of the filter. This is done by first 
placing the zeros on the unit circle at the appropriate 
frequency. Usually the zeros are desired on the unit circle 
to afford as nearly as perfect a stop band as possible. The 
poles are placed inside the unit circle at the proper 
frequency and their distance from the center determines the 
sharpness (quality factor) of the peaks (or pass bands) . 
Figure 3.3a below shows a desired low pass frequency 
response. However, when some insight concerning the pole and 
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zero location is desired, the procedure is as follows: 

1. Place the zeros at the proper position on the 
unit circle for either the high or low pass 
filter characteristic. 

2. Locate poles inside the unit circle at the 
frequency at which they occur on the frequency 
response . 





b. Estimating the Pole and 
Zero Locations . 



FIGURE 3.3 Graphically Estimating Pole 7-ero Locations 
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This procedure can be performed in reverse also. 
For example, if the poles-zero locations are known, the 
frequency response can be estimated by moving 
counter-clockwise around the circle starting at Z= 1+jO. 
The amplitude is increased as a pole is approached and 
decreased as a zero is approached. Midway between the poles 
and zeros, an arbitrary median amplitude value is 
maintained. If the pole is closer to the perimeter of the 
circle the peaks are sharp (high Q) and if they are closer 
to the center of the circle the peaks are not sharp but 
wider (low Q) . 



4 . Integrator a nd Canceller Recursive C omb Fi lter s 



The Recursive Comb filter can be configured as 
either an integrator or a canceller type of filter. The 
difference between these types of filters is simply in the 
relationship between their stop band and pass band widths. 
The integrator is characterized as having a very narrow pass 
band (thus implying a very high Q circuit) and a very wide 
stop band. The canceller is the opposite of this with a 
very wide pass band and a narrow stop band. These 
characteristic responses are shown in figure 3.4. 



33 





Integrator 



Figure 3.4 Recursive Comb Filter Types 
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5. 



Coefficien t De termin atio n 



One portion of this thesis deals with Integrator 

Recursive Comb Filters, and the demonstration of their use 

as a device to improve signal-to-noise ratio of a pulse 

train signal added to white Gaussian noise. The Z transfer 

function is given by equation 1 (given in Chapter 2) and the 

numerator coefficients are known (1, 2, 1 respectively). 

The design of the filter as discussed earlier is a problem 

of determining the denominator coefficients b and b . This 

1 2 

is not a trivial problem since there are very many possible 

combinations of b and b each between ±1 and ±2. Where 

1 2 

then does one begin searching for the proper coefficients 

that will give an integrator Comb filter response and also 
have a very narrow pass band? Due to the massive amount of 
calculations required to plot a frequency response of 
equation 1 the problem lended itself well to to a computer 
iteration. The fortran program shown in Appendix B is a 
programing of equation 1. 

To obtain the frequency response from equation 1, 
Z=Cos(wT| +jSin(wT) was substituted into the transfer 

function resulting in H(jw). This was then programed and 
the only inputs that were required were the numerator and 
denominator coefficients. The output is a plot of the 
frequency response on a linear graph and the amplitude is 
normalited to the maximum amplitude encountered throughout 
the frequency range of interest. This program will 
automatically plot the results on the textronix terminal of 
the CP/CMS time sharing system. 

As noted earlier equation 1 is the transfer 
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characteristic of a digital filter design technique. This 
equation applies equally well to the Sampled-Analog filter 
design technique, although a slight modification must be 
made to account for the sample and hold circuitry which is 
used in the filter to reproduce the filtered signal at the 
output. In, the frequency domain the sample and hold 
characteristic is a SinX/X function which affects the filter 
response in a multiplicative manner. This characteristic 
has been included in the fortran program. One other 
characteristic which was necessary to model the 
Sampled- analog filter was the adjustment of the frequency of 
the sin and cos arguements to give the proper number of comb 
teeth, corresponding to the number of delays which were used 
in each stage of the filter. A second order filter was 
designed with 96 delays per stage. This resulted in a comb 
response with 96 teeth between D. C. and the sampling 
frequency. The 96 teeth are a result of the repetitive 
nature of the bilinear design technique. This 
characteristic was also modeled into the fortran program 
which is shown in Appendix B. 



6 . Fil ter D esig n 

The computer program of Appendix B will give the 
frequency response of the desired Z transfer function, but 
as previously mentioned the coefficients are required input 
data. Thus, a method is needed to determine where to start 
searching for the proper coefficients. This problem was 

[ 6 ] 

answered by the result of a thesis by Stacy V. Holmes 

In his research he derived a means to analize Sampled-Analog 

filters based on the coefficients b and b . In his 

1 2 
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dissertation he introduced a b-plane triangle of stability 



(shown here in figure 3.5 ) . From this graph a rough 

estimation of the denominator coefficients can be obtained. 

although the coefficients which are obtained in this manner 
are very rough figures it does provide a starting point to 
use the previously discussed computer model. Several "b" 
coefficients can be used in the program until the desired 
frequency response is obtained, prior to building the 
filter. Since this study is of recursive Comb Integrators 
and Cancellers, from the graph the b^ coefficients are 

limited to values of approximatly +0.2 to -2.0 where as the 

b coefficient can very from +1.0 to -1.0. This technique 
2 

is , needless to say, time consuming and tedious; however, 

the b-plane triangle of stability decreases the coefficient 
search immensely. 
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Figure 3. 5^^^ b-Plane Triangle of Stability 
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B. COMPOTER AIDED FILTER DESIGN 



One of the intentions in this study is to design a 
Recursive Integrator Comb filter with a very narrow passband 
and a broad stop band and demonstrate its use as a device 
for improving Signal-to-Noise ratio of a pulse train input 
signal (a topic of a later chapter). Suffice it to say that 
the type of filter design desired can be found with 
coefficients that are located in the region below the .05f 

s 

isocontour on the b-plane triangle of" stability shown in 

[ 6 ] . , . . 

figure 3.5 . Several values of coefficients were tried in 

the computer program, these coefficients are as shown in the 
table below: 
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U^lpr. 










Figures 3.6 through 3.25 show the computer results for 
the coefficients listei above. These plots give an 
indication of the results which can be expected 
experimentally from a filter with these coefficients and 
this number of delays. 
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Normalized Amplitude 




Figure 3.6 Comb Frequency Spectrum (Theoretical) 



41 



Normalized Amplitude 







Figure 3.7 Comb Frequency' Spectrum CTKeoretical) 
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Normalized Amplitude 




Figure 3.8 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.9 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.10 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 



b^“*2 "b2~”*d 




Figure 3.11 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.12 Comb Frequency Spectrum (Theoretical) 
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Normalized Amolitude 




Figure 3.13 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.14 Comb Frequency Spectrum (Theoretical) 



49 



Normalized Amolitude 




Figure 3.15 Comb Frequency Spectrum 



(Theoretical) 



50 



Normalized Amplitude 







Figure 3.16 Comb Frequency Spectrum (Theoretical) 



\ 
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Normalized Amplitude 



2 




Figure 3.17 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 



3.18 Comb Frequency Spectrum 



(Theoretical) 
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Normalized Amplitude 




Figure 3,19 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.20 Comb Frequency Spectrum (Theoretical) 
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Normalized Amolitude 




Figure 3.21 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.22 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 
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Figure 3.23 Comb Frequency Spectrum (Theoretical) 
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Normalized Amplitude 




Figure 3.24 Comb Frequency Spectrum (Theoretical) 
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Normalized Amolitude 




Figure 3.25 Comb Frequency Spectrum (Theoretical) 
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IV. 



experimental filter stddy 



A. CIRCUIT CONFIGURATION 



The Recursive Comb filter was connected as shown in 
Figure 4.1 below. 




Figure 4.1 Filter Circuit Configuration 

This circuit configuration is the conical circuit 
realization; using the Bilinear Z transform of either the 
high or low pass S-Plane analog filter transfer function. 
The sample and hold circuit was necessary to give better 
nulls in the frequency response characteristics of the 
filter. If this circuit had not been used, the sumation of 
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the signals at the operational amplifier summer would be the 
addition of a strict analog signal and a sampled analog 
signal. This introduces errors and rapidly deteriorates the 
nulls of the frequency response. With the sample and" hold 
circuit, two similar signals are added in the operational 
amplifier. To further reduce errors, the clock frequency of 
the sample and hold circuit was adjusted to be the same as 
the sampling frequency of the delay elements. The 
coefficients were composed of* carbon resistors with five 
percent tolerences. The resistors were arranged as voltage 
divider networks where the values were adjusted by combining 
resistors in series and measuring the tap voltage (which was 
a portion of the entire voltage across the series resistance 
combinations) . The tap voltage was determined by 
multiplying the total voltage on the entire series 
combinations by the desired coefficient value; and after 
this value was measured at one of the tops, the entire 
combination was placed into the circuit. It was necessary 
to obtain the coefficients in this manner, instead of using 
variable potentiometers for coefficients, since these latter 
devices were found to be frequency dependent (as determined 



in earlier studies 



[ 6 , 7 ]^ 



The delay device was the Reticon SAD-100 serial analog 
delay. In the circuit diagram of Figure 4.1, the total 

delay is represented as 3 . Since the Reticon SAD-100 has 

ninety -six bits of delay, the value of n in this circuit is 
ninety-six clock periods for each delay stage. 



B. EXPERIHENTAL DATA 



The circuit which has been described above was studied 
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to determine its frequeacy response characteristics and the 
dependence of these responses on the denominator 
coefficients. The numerator coefficients were of no 
interest here since it was determined in Chapter Two that 
their relationship must be fixed to give the comb response; 
and, that changing the sign of a^ defines either a high or 

low pass filter (+ for low pass and - for high pass filter) 

characteristic. Therefore, for all data presented here, a 
low pass filter was used. 

The data plan is the same as the plan used for the 
computer design section (table 3) . This was done to compare 
the results of the computer program. Table 3 has been 
repeated here for convenience. 




The experimental data is shown on the following pages in 
Figures 4.2 through 4.17. The two sets of curves for b^ = 
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±1 were not taken since the device was unstable with these 



coefficients. Accompanying each frequency plot is a plot of 
the poles and zeros of the transfer function for the 
coefficients used in the respective response. This was 
included to obtain a better understanding of the reason for 
the shape of the filter teeth. The Z plane plot starts at 
D.C. at Z = 1 + jo and one counterclockwise revolution, 
around the circle, increases the frequency to fs/96 , or one 
complete cycle of the frequency response that it 
accompanies. The clock frequency for the following figures 
was 400 K Hz. 
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Normalized Amplitude 




Figure 4.2 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 




Frequency ( K Hz) 




Figure 4.3 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Figure 4.4 Comb Frequency Spectrum (Experimental) 
and Pole Location 



67 



Normalized Amplitude 





Figure 4.5 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 




Frequency (K Hz) 




Figure 4.6 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Figure 4.7 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 




Frequency fK Hz) 




Figure 4.8 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Figure 4.9 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Figure 4.10 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 




Frequency (K III) 




Figure 4.11 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Figure 4.12 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Figure 4.13 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Frequency (K Hz) 




Figure 4.14 Comb Frequency Spectrum (Experimental) 



and Pole Location 
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Normalized Amplitude 





Figure 4.15 Comb Frequency Spectrum (Experimental) 
and Pole Location 
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Normalized Amplitude 





Figure 4.16 Comb Frequency Spectrum (Experimental) 
and Pole Location 



79 



Normalized Amplitude 





Figure 4.17 Comb Frequency Spectrum (Experimental) 
and Pole Location 



80 



1 . Data Anal y ses 

It appears , from the data in the previous figures, 
that changing the sign of causes the filter response to 

change from a canceller to an integrator type filter. This 

is not completely true. The reason that the filter is 

changing back and forth between canceller and integrator can 

be explained by referring to Figure 3.5 of Chapter 3. In 

this figure, there is an isocontour denuraerated .25fs which 

distinguishes between the two types of filters under concern 

here. The fact that b is held numerically close to this 

1 

isocontour causes the filter characteristic to jump back and 

forth across this isocontour. However, when b =0.2 and b 

1 2 

is greater than 0.6 in absolute value, the filter response 
remains as an integrator since these coefficients fall below 
(integrator side) the .25fs isocontour. There is a large 
difference in the shape of the frequency response teeth just 
described. The changing of the sign of b^ from minus to 

plus causes the poles to change from being real (minus 

coefficients) to being imaginary (plus b^ coefficients) even 

though both curves are integrator responses. 

Another point which was noted was, that as the value 

of b decreases negatively, the slopes of the curve skirts 
2 

become steeper. This was expected, since, as the b-plane 
triangle of stability indicates, when b^ becomes more 
negative, (holding b constant) the value of the isocontour 
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ratio becomes smaller. 



This indicates that the three 



decibel cutoff frequency is decreasing (f was held constant 

s 

throughout the experiment) . This, in effect, causes the 

pass band to decrease, bringing about a sharper response 
curve. This can be seen by comparing Figures 4.3, 4.4, 

4.6, and 4.8 and their counterparts which are plotted with 

b plus. 



Figures 4.14 and 4.16 show the nulls beginning to 
deteriorate quite rapidly. This is due to the shifting of 
the zeros around the unit circle from the Z = ±1 + jo point. 
Since there is a real pole relatively close to the perimeter 
of the unit circle in these two figures, the pole 
immediately becomes effective; allowing some of tbe signals 
at these frequencies to pass through. However, as the 
frequency was increased still further, the second zero is 
encountered and a second null obtained. The deterioration 
became increasingly apparent as the frequency was increased. 
An example of this deterioration is shown in Figure 4.18. 
This figure shows the entire spectrum of the filter response 
from DC to fs = 400 K herrz. Also shown in Figure 4.18 is 
the deterioration of the nulls due to the shifting of the 
zeros and some inefficiency of the Reticon filter. 

The shifting of the zeros may have been due to stray 
capacitances caused by the connecting leads of the circuit. 
This would cause the coefficients of both the transfer 
function numerator and denominator to vary and, in the case 
of the numerator coefficients, cause the zero locations to 
shift off the real axis. This results in an imaginary term 
in the value of Z, causing a shift around the unit circle. 
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Figure 4.18 Comb Frequency Spectrum (0-400 KHz) 



c. 



FILTER SELECTION 



As will be discussed in the following chapter, an 
integrator recursive comb filter is desired with very narrow 
pass band teeth to demonstrate the Signal-to-Noise ratio 
improvement ability of the comb filter. The type of 

response chosen was the filter with b = -0.2 and b = -0.6. 

1 2 

This filter has a very narrow three decibel pass band 
(approximately 417 hertz wide) and a wide stop band 
(approximately 3.91 K hertz wide). This is shown in Figure 
4.6. Figure 4.19 shows the entire frequency response from 
DC to 400 K Hertz. The zeros are relatively unaffected at 
the higher frequencies although the amplitudes of the teeth 
are increasing between 80 K hertz and 110 K hertz. This was 
possibly due to the frequency dependence of the b 
coefficients since, as explained earlier, these coefficients 
affect the teeth shape. However, for the purposes of the 
next chapter, this was acceptable. The discrepancy in the 
teeth amplitude at approximately 200 K hertz is due to 
instrumentation errors only. 
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Figure 4.19 Comb Frequency Spectrum (0-400 KHz) 



V. 



APPLIC A TI ON S OF THE INTEGRATOR TYPE RECORSI_VE COMB 
FILTERS U SI3MAL;TO-NgiSE RATIO IMPROVEMENT 



A. INTRODOCTION 

Signal-to-Noise ratio (SNR) iniprovement is a aajor 
concern in radar and communication design problems. In this 
section an attempt has been made to use a Charge Transfer 
Device (CTD) to improve the SNR. The input signal that was 
chosen was a pulse train. A pulse train was chosen since it 
is commonly used in radar applications as well as the fact 
that it has a discrete frequency spectrum with many 
frequency components. This latter fact demonstrates the 
intent of this section more explicitly, as will be seen 
later. 

B. FOURIER ANALYSIS AND FILTER CONSIDERATIONS 

The result of a Fourier analysis of a pulse train is 
given by equation (1) below and its development is shown in 
Appendix A. 




-jkw_(to + t/2) SinkWot/2 

) ( 

kwo t/2 



1 



where 



t =pulse width 
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r= pulse repition frequency (PRF) 
V= pulse amplitude 
w= 2PI/T 



The line spectrum is shown in Figure 5.1. Note that the 

amplitude is a function of the familiar SinX/X curve (where 

X = kwt/2 for this explanation ) and that all the even and 

odd harmonics are present since the pulse train in this 

situation is neither an even nor an odd function of time. 

Further note that the value of C is directly related to the 

k 

pulse width (t) and inversely proportional to the pulse 
period T. If the pulse period (T) increases, the band width 
of the spectrum ( band width is from 0 hertz to 1/T) does 
not change, but the number of spectral lines between freg.= 
0 and freq.= 1/T increases. If this pulse train were 
transmitted through a very noisy environment it would be 
most difficult, if not impossible, to detect it unless 
special filtering techniques were employed. However, if the 
frequency spectrum of the pulse train could be matched with 
the frequency spectrum of a filter which will attenuate all 
other frequency components in between the desired signal 
components, then the system noise will have minimum effect 
on the transmitted/received signal. The problem then is to 
find a filter which has the same frequency spectrum 
characteristics as that of Figure 5.1. A filter with this 
type of spectrum is the Integrator Recursive Comb filter. 
This filter is periodic in the frequency domain and is 
capable of nearly perfect attenuation within the stop band. 
Since the spectrum of the pulse train has discrete frequency 
components over a very wide frequency range a very narrow 
comb filter pass band response is desired to capture the 
signal's spectral components and at the same time attenuate 
as much of the undesired spectral components as possible. 
The characteristic response of the Integrator Comb filter is 
such a response. The characteristic response of the 
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Integrator 
response can 
techniques, 
since there 
conversion 
comb filter 
amplitude o 
frequency f 



Comb filter is shown in Figure 5.2. This 
be produced very easily using sampled analog 
Sampled analog techniques are highly desirable 
is no analog- to-digital or digital- to- analog 
required. Also, the frequency spectrum of the 
can be shifted without changing the relative 
f the pass band, simply by changing the sampling 
. Since the comb filter frequency response has 



s 

a finite band width pec tooth, some of the noise spectral 



components will pass through the device with the signal's 
spectral components. The major portion of the noise should 
be attenuated since the pass band widths are much narrower 
than the stop band widths of the filter response. 
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Normalized Amplitude 
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Figure 5.1 Pulse Train Frequency Spectrum 
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Figure 5.2 Integrator Comb Filter' Response 
C. DEVICE SELECTED 

The device selected to compose the filter was the 
Reticon SAD-100. This device is a monolithic integrated 
circuit which provides storage of analog information. The 
device has 96 bits of delay and provides a first-in, 
first-out storage for these 96 samples. The specification 
sheet on the SAD-100 indicates that there are 100 bits of 
delay, two of which are used for input and output junctions 
leaving 98 storage cells. Earlier work with this device had 
indicated that there were 96 storage cells. The data which 
was obtained throughout this study proved this to be 
accurate. 

Some of the important characteristics of the SAD-100 are 
given below; 



90 



Maximum Delay at 4 KHz sampling rate 25 msec 

Maximum Sampling Hate 12 MHz 

Dynamic Range 70 db 

Minimum Clock Frequency 2.0 KHz 

Storage Cell Capacitance V. 2 pf 

Input Capacitance 8.0 pf 

Output Capacitance 8.0 pf 



The coefficients which were used are: 



a =1.0 
0 

a = 2 . 0 b 

1 1 

a = 1.0 b 

2 2 



- 0.2 



- 0.6 



These coefficients give the most desirable frequency 
response as shown in Figure 5.3. The stop bands are nearly 
perfect zeros and the pass bands very narrow with very sharp 
slopes on both sides of their respective center frequencies. 



The filter was operated at a sampling frequency of 390 K 
hertz during the measurements. Thus, if the pulse train is 
to pass through the filter unattenuated and without 
distortion, the spectral components of the pulse train must 
fall directly on the pass bands of the filter response. 
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Figure 5.3 Comb Filter Frequency Spectrum (Experimental) 



The pass bands of the Comb filter are periodic with 

respect . to the sampling frequency (f ) divided by the 

s 

number of delays in the SAD - 100. The filter's first tooth 
center frequency is given by equatioa 2 below: 

f = f /Number of Delays 2 

1 s 

This calculation gives a center frequency of the first 

tooth of the filter at f =4063 hertz. Since the following 

1 

teeth are repetitive with respect to this frequency, the 
following teeth center frequencies are given by M*f , where 
H is an integer equal to 1, 2, 3, ... 



For the signal to pass through the filter, the spectral 
components of the signal and those of the filter must 
coincide. For the pulse train to coencide with the filter 
frequency spectrum, it must have a pulse repitition 
frequency of 4063 Hertz or some integer multiple of this 
frequency. The frequency selected was 4063 Hertz and its 
frequency spectrum is shown in Figure 5.4, Figure 5.5, and 
Figure 5.6. Figure 5.5 and 5.6 are an expanded scale of the 
input signal , whereas figure 5.6 is the entire spectrum 
from DC to 450 KHz. A comparison of figures 5.3 and 5.4 
will reveal that the spectrums are matched perfectly. The 
time domain representation of both the input and output 
pulses are shown in Figure 5.7. In this figure only one 
input pulse is shown, and above it, one output pulse. The 
slight distortion in the output pulse is due to the loss of 
some of the higher frequencies resulting from an imperfect 
frequency match. This loss is caused by slight drifting of 
the pulse generator from the selected frequency, as well as 
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the zeros of the stop band at the higher frequencies rising 
above zero. Figure 5.8 shows a series of input and output 
pulses. Figure 5.9 demonstrates the signal attenuation 
when the input signal frequency spectrum is not matched to 
the filter spectrum. The input frequency in this case is a 
3.56 K hertz pulse train ( shown on the bottom of figure 
5.9); note that the output signal is completly Lost (shown 
in the top of figure 5.9). 



; 
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Figure 5.4 Pulse Train Frequency Spectrum (PRF=4063 Hz) 







Figure 5.5 Input Pulse Train Frequency 
Spectrum (Experimental) 




Figure 5.6 Input Pulse Train Frequency 
Spectrum (Experimental) 
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Figure 5.7 Input Pulse 




Figure 5.8 Filter Input and Output 
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Figure 5.9 Cancellation of the Input Signal 
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The input noise which was used was produced by General 
Radio Company and was a type 1390-5 random noise generator. 
This device produces a white Gaussian noise over a frequency 
range of up to 5 M hertz. The frequency spectrum of the 
Gaussian noise was plotted on the X - I recorder and is 
shown in Figure 5.10. This graph is plotted from zero hertz 
to 500K hertz. The slight dip from approximatly 100K hertz 
to 300K hertz is due to the noise generator characteristics. 

The noise as seen in the frequency response of Figure 
5.10 was passed through the Reticon integrator Comb filter. 
As to be expected, the filter frequency response should be 
similar to that of figure 5.3. The actual noise output 
frequency response of the filter is shown in Figure 5.11. 
This plot is from zero hertz to 50K hertz. A comparison of 
this and figure 2's response reveals that the bass band has 
widened and that the peaks have become jagged. This is due 
to the noise having a very low input amplitude of 0.6 volts 
RMS. Thus, the product of the very low input signal and the 
filter pass band tooth on a normalized basis causes the 
small amplitude skirt portion of the pass band to remain 
approximatly the same as it was before but; the product of 
the very large center portion of the filter pass band with a 
very small input spectral components causes the large center 
frequency amplitude to decrease considerably giving the 
effect of widening the pass band of the filter on a 
normalized basis. However, this in effect is not the case. 
(Note that all of the actual curves plotted have been 
normalized to ease the comparisons which have been made.) 
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Figure 5.10 Input Noise Spectrum 




9pnq.T'[dutv poziiBuiaoM 



101 



Figure 5.11 Output Noise Spectrum 



Figurs 5.12 graphically illustrates the attenuation of 
the noise signal with the input noise spectrum and output 
noise spectrum plotted on the same graph. The graph is 
plotted from zero to 50K hertz. The noise input spectrum is 
the top curve and the output spectrum is below it. This 
same input and output noise are shown in Figure 5.13 in the 
time domain. Note that although much of the noise is passed 
through the filter, much of it has been cancelled out. In 
Figure 5.13, the bottom is the input noise and the top, the 
output noise. 

/ 
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Figure 5.12 Comparison of Input and Output Noise Spectrum 




Figure 5.13 Time Domain .o£ Input and Output Noise 
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D. 



SIGNAL-TO-NOISE RATIO MEASOREMENTS 



It has been determined at this point that the Integrator 
Comb filter will act well in cancelling much of White 
Gaussian noise. But, can a signal be detected well enough 
to determine its presence ; and, if so, what type of signal- 
to -noise ratio improvment can be obtained? These are very 
important questions and will be pursued next. 

The frequency spectrum of the input pulse train and 
noise was added with a resistive network as shown in Figure 
5.14 below. 



1000 ohms 




Figure 5.14 Network for adding Signal and Noise 



T05 



The resistances were necessary to prevent the noise from 
feeding into the signal source as well as to prevent the 
signal from feeding into the noise source. The l.OK ohm 
resistor to ground was necessary to prevent the relatively 
low input impedance of the filter from changing too 
drastically . 

Figure 5.15, Figure 5.16, Figure 5.17, and figure 5.18 
are photographs of the signal plus noise input and output. 
The input signal was decreased from 2 volts peak to peak 
(figure 5.15) to 0.1 volts peak to peak ( figure 5.18). 
This was done to give an indication of how well the filter 
will reproduce the signal even though it is buried in a 
relatively high noise environment. As seen in figure 5.18 
the output pulses ( bottom trace ) are much more noticable 
than the input pulses. 

Figure 5.19 demonstrates the cancellation of the input 
signal, which is added to the Gaussian noise, when the input 
frequency is not matched to that of the filter. Note that 
the output trace (bottom ) is only composed of filtered 
noise whereas the top trace ( input ) is composed of both 
the signal and noise. 

The previous data gives a visual indication that the 
signal is reproduced and also gives some indication of how 
well the signal is reproduced. However, a more precise 
method in answering the question of how well the signal is 
reproduced is to compare the signal's input and output 
frequency spectrum's. The following photographs show the 
frequency spectrum comparisons of the input and output. 

Figure 5.20 is the frequency spectrum plot of the 
output pulse train alone. It covers a frequency range from 
zero to 450K hertz. Note that there is a large impulse at 
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390K hertz. This is dae to the sampling clock PRF feeding 
through the filter; not a characteristic of the pulse train. 
Further note that the spectrum decreases to approximatly -67 
dB at its minimum value which occurs at 390. OK hertz. This 
is as predicted since the Fourier analysis (as seen in 
appendix A ) predicts that the spectrum will have a SinX/X 
envelope with the first null occuring at the reciprocal of 
the pulse width ( i.e. 390. OK hertz ). Figure 5.21 shows an 
overlap of both the input and output spectrums. The two 
curves are very much alike with the exception that the 
output spectrum ( top curve ) is approximatly 10 dB higher 
than the input spectrum (bottom curve ) . This is due to the 
filter having a gain of approximatly 10 dB. Figure 5.22 is 
the same output spectrum of figure 5.20 except that the plot 
extends only to approximately 20K hertz. Finally, Figure 
5.23 shows the input and output plots overlapping. Again 
note the increase in the peak amplitudes due to the gain of 
the filter. Thus this comparison of both the time domain 
and frequency domain proves that the pulse train is 
reproduced exactly. 



/ 
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Figure 5.15 Signal Plus Noise 




Figure 5.16 Signal Plus Noise 
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Figure 5.17 Signal Plus Noise 




Figure 5.18 Signal Plus Noise 
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Figure 5.19 Cancellation of Signal Plus Noise 




Figure 5.20 Output Signal Frequency Spectrum 
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Figure 5.21 Comparison of Input and Output Spectrums 




Figure 5.22 Output Frequency Spectrum 






Figure 5.23 Output and Input Frequency Response Comparison 
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E. FREQUENCY CONSIDERATION" 



To this point there has been no comparison of signal 
plus noise in the frequency domain. The following data will 
show that the input signal is uneffected by a noisy 
environment and is reproduced exactly by the filter. 

Figure 5.24 is a frequency spectrum of the input signal 
plus noise from zero hertz to 50 K hertz. The spectrum 
shown is a continious one since Gaussian noise has a 
continious spectrum; but , in this graph the pulse train 
spectral components are added to the noise spectrum causing 
the spectral lines to protrude higher than the noise 
spectral amplitudes. For a signal - to - noise ratio 
improvement it will be necessary for the Integrator Comb 
filter to attenuate the noise spectral components. This is 
in fact what has been done as seen in Figure 5.25. All the 
noise spectral components between the signal spectral lines 
have been attenuated to nearly zero. Note, however, that 
there is still noise passed through the filter since the 
pass band teeth have some finite width. Also, the width of 
the output spectrum teeth appear to be much wider than the 
filter pass band teeth. This is due, as explained earlier, 
to the small signal levels. Thus from the previous data, 
the signal is in fact reproduced exactly since all the 
spectral components of the pulse train signal are passed 
through the filter but also along with this reproduced 
signal there is noise passed through; however, the noise has 
been reduced considerably. 
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Figure 5.24 Pulse Train Plus Noise Frequency Spectrum 
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1-igure 5.25 Filtered Pulse Train Plus Noise 



F. 



SNR IMPROVEMENT USI N3 A SECOND ORDER FILTER 



It is evident at this point that there has been an 
improvement in signal - to - noise ratio, since the signal 
has been reproduced exactly and much of the noise has been 
attenuated. The question still remains as to how much 
improvement has been obtained. The answer to this question 
was obtained by making voltage measurements with an RMS 
voltmeter. The measurements were as follows: 

Pulse Train Voltage (RMS ) .104 .880 

Noise Voltage ( RMS ) .665 1.15 

Filter Noise Voltage < RMS ) - .033 

Since these are voltage measurements the signal-to-noise 
ratio is calculated using equation 3 below. 

SNR(dB) = 20x(Log V - Log V . ) 3 

10 sig 10 noise 

Thus the values are as given below: 

SNR = - 16. 12dB 
in 

SNR = -2.32 dB 
out 

The difference of the output SNR and the input SNR is 
the improvement. For these conditions it is; 

SNR = 13.79 dB 

improvement 



Although a 13.79 d3 SNR improvement is not a great feat 
the fact remains that, for a first generation application of 
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the Reticon SAD - 100 as an Integrator Recursive Comb 
filter, it has performed well. Room for improvement in its 
operating techniques as well as the filter design techniques 
still remains. 



G. FIRST ORDER FILTER SNR IMPROVEaENT CORVES 



The work which has been done above was all with a second 

order integrator filter. The following data was taken with 

the Reticon SAD-100 operating as a first order filter having 

numerator coefficients as: a = a = 1.0. The b 

0 1 

coefficient was variable for this demonstration. The 



transfer function for this filter configuration 
below as equation 4. 



is shown 



H(z) 



1 + Z -1 
1 + b, Z-i 



4 



The purpose of this section 
of the SNR improvement ability of 
The experimental data is shown in 
plotted and is presented below in 



was to obtain some insight 
the first order filter. 
Appendix C. This data was 
figure 5.26. 
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rigure 5.26 SNR Improvement vs b. Coefficients 



The graph of figure 5.26 above is plotted with the b 

1 

coefficient on the abscisa as the independent variable and 

the SNR improvement as the dependent variable. the data was 
plotted in this manner since the b^ coefficient determines 

the shape of the filter passband teeth as well as whether 

the filter is a canceller {+ b’s) or an integrator (- b's) . 
Thus by varyino b^ the filter will either pass more or less 

noise component through with the pulse train spectral 

components giving a variation in the SNR improvement. 

Three curves are plotted each with varying SNR 

input 

values. Although curves A and B only differ by 0.4 db, the 

input signal and noise voltages doubled in curve B giving 
approximatly the same input SNR. This was done to see the 
effects the input anlplitude had on the amount of 
improvement. Curve C is a plot with the input SNR 
approximatly three times the value of either curves A or 3. 
It was interesting to note that all the curves are within 
2db or less of each other as seen in figure 5.27. This was 
an expected result since the filter pass band teeth are the 
same for the respective b coefficients regardless of the 
input signal amplitude. The SNR improvement is much higher 
for negative b coefficients because of the filter teeth 
becoming much narrower (i.e. a higher Q integrator circuit) 
and attenuating more and more of the noise. Since the three 
curves of figure 5.26 are all relatively close to each 
other, it was felt that the variations are due to 
instrumentation measurement errors and the fact that the 
noise is such a highly random signal that the measurement of 
the noise changes slightly each time a new measurement is 
made on the same signal. These small changes in 



measurements cause a noticable logarithmic variation due to 
the nature of the logarithmic conversion characteristic. 




Figure 5.27 Band o£ Data Variation 
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VI. APPLICATI ON OF THE IN TEG RATOR TYPE REC URSIVE COilB 

2i SE LE CT ION OF PULSE TRA IN 



A. RADAR 



Now that the Integrator Recursive Comb filter has been 
demonstrated as a device with potential to improve 
Signal-to-Noise ratio, to what practical use could such a 
device be put? One application which has been discussed in 
the literaxure (Recent Advances in the Synthesis of Comb 
filters) is their use for Radar application for video 
integrators. They are used in this application to 
accomplish linear addition (integration) of the radar 
signals which are consistant from sweep to sweep and at the 
same time the random signals ( noise ) will add in a root 
sum -square form. The net result is an improvement in the 
SNR as shown in the previous section. 



3. ELECTRONIC COUNTERMEASURES 



Another application which is proposed by this report is 
the use of the Integrator Recursive Comb filter as a Pulse 
Sorter of Radar pulses to speed up the acquisition of Radar 
signals in an Electronic Countermeasures situation. This 
concept of a pulse Sorter has been demonstrated and the 
results are shown in Figure 6.1. The systems pulse trains 
added together make up the first four traces. There are four 
separate traces to show each pulse train as a stationary 
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• I 



4 







signal on the scope. This was done by synchronizing the 

oscilloscope separately with each input signal and 

reexposing the same negative to the scope trace, after 

synchronizing and shifting the position of the stationary 

pulse train. After the first exposure, the filter output 

was moved off the scope screen, the next pulse train was 

synchronized to the scope, (to keep it stationary) the 

entire trace was moved vertically down on the screen, and 

the film was reexposed to the screen. Thus, the same 

negative was exposed four times. The frequencies of each 

pulse train from top to bottom are: f = 8.46K hertz, f = 

1 2 

1000 hertz, f = 2000 hertz, and f = 1500 hertz. The last 
3 4 

trace is the filter output of the captured pulse train (f^ 

in this situation) . Note the correspondence between the 

pulse locations in the top trace f^ and the output trace 

(bottom trace). This is the captured signal. The others do 

not pass through the filter since they fall in the filter's 
stop bands and are attenuated. Thus, by using the Integrator 
Recursive Comb filter, a particular pulse was sorted from 
among several and interference reduced from non-synchronous 
pulse train. 
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Figure 6.1 Pulse Sorting Illustration 
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• I nterfere nce considera tions f £om non- synch ro no us 
trains 

Ideally, the Eeticon filtering will not pass through 
the system at all. But, there are some frequency 

relationships between the desired passed pulse train and the 
undesired signal so that the filter will allow for passage 
of some of the undesired components. The components that 
will pass are those of the signals that have a whole 
multiple spectral component of the filter's clock frequency. 
For example, the frequencies which were used above (f^ = 

8.46 K hertz, f = 1000 hertz, f = 2000 hertz, and f = 1.5 

2 3 4 

K hertz) , f , f , and f have spectral components that are 
2 3 4 

8f = 8Khertz, 4f = SKhertz, and 8f = 12Khertz and will 
2 3 4 

fall within the pass band of the filter teeth. The effects 

of these components will cause some (as yet undetermined) 
distortion. The distortion was not noticable in the data 
taken in this thesis. An indepth study of this distortion 
would be most interesting as a thesis chapter in the future. 
Figure 6.2 below shows the spectral components that fall 
within the filter pass bands and very close to the captured 
signal ( f ^) • 
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Signal Frequency = ]000 Hz 






Freq. (KHz) 
Signal Frequency = 8460 Hz 



4.23 8.46 ]2.69 ]6.92 

Freq. (KHz) 

Figure 6.2 Interferring Spectral Components 
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Frequenc y S 2 ectrjim Considerations 

The previous figures have demonstrated that the 
filter can select a specific pulse train signal and pass it 
while attenuating the others. The signal which is selected 
by the filter is reproduced at the output; but, how 
accurately is the signal reproduced? Figure 6.3a shows, on 
a linear plot, the line spectrum of the input pulse train. 
This is the signal that was eventually passed through the 
filter to give the output trace in Figure 8.1 above. 
Although the filter will pass multiples of signals starting 
at 406K /96 Hertz, the primary frequency selected for this 
demonstration was 8.46 K hertz. With this frequency, every 
other tooth of the comb filter frequency response will 
contain a spectral component of the 8.46 K Hertz input 
signal. This was done to ease the comparisons which are 
made later in this report in determining the filter's 
selectivity. Figure 6.3b shows the output of the filter 
using the input signal as shown in Figure 6.3a. Note that 
there are some very small spectral components at 
approximately 4K Hertz, 12K Hertz, 20K Hertz, etc. These 
components are that portion of the comb filter's teeth pass 
bands that are not filled with spectral lines of the signal. 
Instead, there are some noise components (possibly filter 
generated) which fall into these unused teeth and are being 
passed through. In comparison of these two signals. Figure 
6.4 shows an overlay of both the input and output line 
spectrums. It is important to note that the signals' 
spectral lines match perfectly and at the proper frequency. 
This shows that there is no frequency distortion as the 
signal passes through the filter, but does not give any 
indication of phase relationship. Note, however, that the 
spectral amplitudes are very different. This amplitude 
change is a result of the filtering system gain. 
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Figure 6.3a Output Frequency Spectrum 




Figure 6.3b Input Frequency Spectrum 
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Figure 6.4 Comparison of Input and Output Spectral Components 
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In summerizing the observations of Figures 6. 1,6. 2, 6. 3, and 
6.4 is noted that; first, the filter reproduces the input 
signal without frequency distortion; second, the filtering 
system amplifies the signal; and third, the filtering system 
introduces some small amount of noise into the signal. Now 
that it is known that the filter can reproduce the input 
pulse train with an added benefit of amplifying it; what 
will be the effects of two input pulse trains with equal 
PRF's and amplitude be on the filter? The following data 
will answer this question. The spectral components of the 
8.46K hertz pulse train are shown in Figure 6.5a. It was 
expected that if a second pulse train of an equal PRF was 
injected at the filter input, the spectral components would 
simply add. This effect was demonstrated and is shown in 
Figure 6.5b. 




Figure 6.5a Spectral Output with Two Input Signals 
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Figure 6.5b Spectral Output with Two Passed Signals 



/ Figure 6.6 shows the overlap of 
Figure 6.5. Note that the two signal 
simply add and that their frequency 
unchanged. 



both the curves of 
spectral components 
relationship is 




Figure 6.6 Comparison of Figures 6.4 and 6.5 



t 
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The effects of hariaonically related signals were 
demonstrated in the above data. These effects were used to 
aid in determining the frequency selectivity of the filter. 
This was done by using two input pulse trains which were 
harmonically related; 4.44K hertz and 8.88K hertz. The 
spectral components of the 8.88K hertz signal are 8.88K 
hertz, 17.76K hertz, 26.64K hertz, etc. and those of the 
4.44K hertz signal are 4.44K hertz, 8.88K hertz, 13.22K 
hertz, 17.76K hertz, etc. Thus, the line spectrum of the 
sum of these two signals will reveal that the spectral 
components ox the 8.S8K hertz signal should be larger than 
those of the 4.44K hertz ones because they are an addition 
of both the 8.88K hertz and 4.44K hertz frequency 
components. Figure 6.7a shows this effect on the linear 
spectrum plot. This spectrum is the addition of two sign-als 
(8.88K hertz and 4.44K hertz) and was used as the input to 
the filter. Figure 6.7b shows the filter output. In both 
pictures, the larger amplitude components 'represent the 
summation of the 4.44K hertz and 8.88K hertz spectral 
components; and the smaller spectral lines between are only 
the 4.44K hertz spectral lines. 
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Figure 6.7a Input Spectral Components 
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Figure 6.7b Output Frequency Spectrum with Two Input Signals 
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The differences in spectral component amplitudes was 
used in determining at what frequency the 8.88K hertz signal 
was attenuated and not present at the filter output. Then, 
by taking the difference of the passed signal frequency and 
the unpassed frequency , the filter selectivity could be 
determined. Figure 6.8 shows a comparison of the output of 
the 4.44K Hertz and 8.78K Hertz input signals. Note that 
the components of the 8.8K Hertz signal are deteriorating 
and decreasing. Decreasing the signal further to 8.59K 
Hertz as shown in Figure 6.9 shows that the original 8.88K 
Hertz is nearly gone with the exception of some 8. OK Hertz 
and 16K Hertz components. The filter output at 3.59K Hertz 
was very distorted and the pulse train was undetectable 
visibly on the oscilloscope. Thus, using this frequency as 
the limit of the filter's passable frequency, results in a 
selectivity of 290 Hertz on either side of the passband 
center frequency. 
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Figure 6.8 Comparison of Two Harmonically Related Signals 




Figure 6.9 Limit of Frequency Selectivity 
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c. 



PROBLEHS ENCOUNTERED IN DATA TAKING: 



One of the major problems encountered during data taking 
was the problem of the sampling frequency clock drifting as 
much as 10K Hertz. Slight adjustment in the circuitry was 
continuously necessary. Also related to the clock problem 
was the problem of keeping the clocks in synchronization 
when the filter was connected in the second order design. 
However, once synchronised they generally remained so. 

The pulse generator was another source of difficulty in 
that the selected frequency for both the pulse Sorter and 
the SNR improvement demonstration would drift enough to 
cause the output to be distorted appreciably. 
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VII. CANCELLER APPL ICAT IONS AND DE MONS TRATION 



A. BACKGROUND 



The most common use of the recursive canceller comb 
filter is as a cancellation circuit for a Moving Target 
Indicator (MTI) radar. In this type of system it is 
im-portant to attenuate the unwanted signals (called clutter) 
from the radar return spectrum. This clutter results from 
stationary targets such as; mountains, buildings from the 
local surroundings, and moving targets such as the leaves on 
trees which are moving with the wind , the trees themselves 
swaying with the wind, ocean waves and swells, etc.. It is 
highly desirable that this narrow band of clutter signals be 
cancelled and the desired signal targets pass through the 
filter and eventually give a target presentation on the 
radar scope. The recursive canceller comb filter is ideal 
for this application. Since it is periodic, a majority of 
the narrow frequency-banded clutter, which is periodic due 
to the pulse repition frequency (PRF) of the radar signals, 
will be cancelled. Also, the canceller frequency spectrum 
characteristic has a wide pass band and will pass the 
doppler shifted moving target. (unless it falls in the 
filters stop band which is then at the radars' blind speed). 
Complete cancellation of the clutter will not occur, 
however, since the stop band of the filter is generally not 
as wide as the clutter spectrum. If the- clutter is a single 
spectral component, (stationary targets) then the filter 
will cancel it. 
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B. CANCELLER FIITFR DEMCNSTRATION 



The technique of cancelling an unwanted frequency and 

passing a wanted signal was demonstrated with the F.eticon 

SAD-100 used previously in this study. The filter was 

connected as a first order comb filter. The clock frequency 

was reduced to its lowest possible value of 104.6 KHz. 

Although the filter is capable of operating at a clock 

frequency as low as 2 KHz, this frequency was not achievable 

without physically changing a capacitor end a variable 

resistor on the printed circuit board that the SAD-100 and 

its associated circuitry was mounted on. The clock 

frequency of the filter was reduced to allow cancellation of 

as low a frequency as* possible to simulate the previously 

discussed clutter spectrum. However, at this clock 

frequency, the lowest frequency that could be- cancelled was 

DC and multiples of 1.09K Hertz. The frequency response of 

the filter at f =104.6 KHz is shown in figure- 7.1. The stop 
^ c 

band center frequencies are at whole number multiples of 
1.09 KHz as seen in figure 7.1. For demonstration purposes 
the unwanted and desired signals were arbitrarily chosen to 
be 2 KHz and 13.1 KHz. Figure 7.2 depicts, in' the top 
oscilloscope trace, the addition of the two signals at the 
filter input. The bottom trace shows the filter output of 
the 13.1 KHz signal. Figure 7.3 shows the passage of 13.1 
KHz signal and a 1.5 KHz signal. This demonstrates that the 
filter will pass two simultanious signals and that if the 
unwanted signal is not at the stop band frequency, it will 
not be attenuated. 
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Normalized^ Amplitude 



fc = 104. 6K bj * 4.4 ^ -1 




Frequency (K IIz) 



Figure 7.1 Frequency Response 



f , , =104.6 

clock 



KHz 
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Figure 7.2 Demonstrating the Cancellation of a Signal 




Figure 7.3 Demonstrating the Cancellation of a Signal 
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VIII 



SUMMARY and RECOMM MDATIONS 



Results of this thesis are summarized ia three areas: 
Experiment, Theory, and Applications. 



A. EXPERIMENTAL RESULTS 



Sampled analog recursive comb filters using a first or 

second order CTD delay line have been studied experimentally 

in detail. The filter circuit used for this study was the 

canonical circuit universally used for the digital recursive 

filters. The feedforward coefficients in this circuit (a , 

0 

a , and a ) are restricted to a mutual relation of 1, ±2, 1 
1 2 

in order to achieve perfect separation (infinite 

attenuation) between neighboring comb teeth. However, 
experimental results revealed that the nulls between comb 
teeth cannot be maintained. Deterioration of these nulls 
results with increasing frequency. To combat this problem, 
it was found that a sample and hold circuit was necessary at 
the input of the filter. The nulls of the filter frequency 
response improved in most cases giving nearly perfect 
attentuation between the neighboring comb teeth. 

The major portion of this thesis dealt with the 
integrator comb filter while a companion thesis (written by 

C 11 ] 

Arif Ejas ) concentrated on the canceller comb features. 
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The Reticon-96 bit delay SAD-100 was the device used to 
implement the comb filter. This device’s characteristics 
were studied extensively and the measurements taken are 
contained in Chapter 4. This study revealed that sampled 
analog recursive comb filters deviated from the ideal comb 
characteristics in that; 

1. The peaks of the comb filter frequency response 
vary with frequency due to the Sin X/X characteristic of the 
sample and hold circuit, and also with the frequency 
dependence of the numerator and denominator coefficients of 
the filter transfer function. 

2. A sample and hold circuit was necessary at the 
front of the canonical filter circuit to improve the nulls 
in the frequency response. 



B. THEORETICAL ANALYSIS 



To model the sampled analog recursive comb filter from 
the existing digital recursive filter theory, it was 
necessary to consider the following: 

1. The feedforward coefficients a , a , and a must 

0 1 2 

be related by the ratio of 1, ±2, 1 respectively. This was 
necessary since the coefficients will place a zero on the 
unit circle in the 2-plane at either Z = +1+j0 or 2 = -1+jO, 
resulting in the desired nulls. 

-1 

2. The single delay 2 used in digital recursive 

-n 

filter theory, must be changed to 2 where n is the number 

of delay stages in the device. This accounts for the n comb 

teeth from DC to the sampling frequency (f ) . 

s 
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3. The filter coefficients are all dependent on 
frequency due to several causes: sampled and hold effect, 
frequency rolloff due to limited bandwidth of electronics, 
and charge transfer inefficiency effect. 

In summary, the theoretical transfer function for a 
second order sampled analog recursive comb filter is: 

a^ (f) ♦ a, (f) Z"'' + aj (f) 

H (z) 

1 + b, (f) Z"^ + bj (f ) Z"^" 

where a = 1, a = ±2, a = 1 at low frequency. 

0 1 2 

Dsing this modified theory, the theoretical analysis and 
experimental measurements are in close agreement. However, 
if proper frequency dependence of the filter coefficients 
can be obtained, the agreement between theory and 

experimental measurements will be improved. 



C. APPLICATIONS 



Three applications have been demonstrated, two for the 
integrator comb filters and one for the canceller type comb 
filters: 

1. Integrator comb filter 

a. Signal to noise ratio improvement 

b. Pulse sorting 

2. Canceller comb filter 

V 

a. Deleting unwanted sinusoids 
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D. CONTINOATION OF THIS RESEARCH 

Suggestions for future work are listed below: 

1. Study the frequency dependence of the 

coefficients. This study would involve more 
computer simulation and matching of actual 
data and the computer plots such as those 
contained in Chapter 3. 

2. Physically implementing the pulse sorter at the 

front end of a wide band receiver to evaluate 
its use under physical "real world" 
conditions. 

3. Implement the Reticon filter with a chirped 

clock frequency and evaluate its performance. 
This type of device would have application 
in a chirped or jittered PRF radar. 

4. Determine the input-output phase relationship 

of the pulse sorter of Chapter 6 to determine 
if there is any phase distortion encountered 
by the filter. 

5. Determine the affects of the unwanted spectral 

components that fall in the filters pass 
bands on the desired pulse train spectral 
components for the pulse sorter discussed 
in Chapter 6. 

6. Theoretically compare actual data with the data 

of figure 5.26 based on area under the comb 
teeth versus the response attenuation area. 
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APPENDIX A 



POLSE TRAIN FOORIER ANALYSIS 




Figure A.l Pulse Train 
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Rewriting 



C = 



k 

C : 



k 



V/- jkWoT( 



- jkw^ (to +t) 
e 



- jkwo to 

e ) 



-jkWo (to/2 + t/2) -jkWo (to/2 + t/2) 

V/-jkWoT( e e 

-jkWo (V2 +t/2) -jkWo(V'2 - t/2) 

- e e ) 



-jkWo (t/2 +t/2) {y2 + t/2) 

V/-]kwoT ( e ) ( e ^ 

-jkWo (V2 ' t/2) 

- e ) 



-U a,/ (V2 +t/2) -jkWoV2, ^ -jkw„t/2 

= V/- 3 kWo T ( e ) ( e ) ( e 

+jkWot/2 

- e ) 



using the Trig. Identity: 



jX -jX 
SinX = (1/2j) ( e - e ) 



-jkWo(to+ t/2) +jk^t/2 -jkWot/2 

= 2V/kH;r( e ® ” ) (e ^ - e ° ) /Ij 



-jkWo (to+ t/2) 

= 2V/kWQ ( e ) ( SinkWot/2) 



X = kw t/2 
0 



Therefor : 
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Solving for the absolute value of C : 

k 



= (Vt/T) ( SinX/X) 



Similarily solving for the DC value of the pulse train 

T 

C = a = 1/T |f(t)dt = 1/T \ Vdt 

0 0 J J 

o to 



C = (1/T) Vt 

0 



to 

to 



Therefore; 



C = (1/T) 7 ( t +t - t ) 
0 0 0 



C = (1/T)Vt 

0 
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APPENDIX B 



COMPlfTER PROGRAM FOR PLOTTING H< Z > 



READ (5,200) AO, A1 , A2, B1 , B2 
200 FORMAT (7F6. 3) 

PI=3. 141 592654 
DEL=0. 1 
HMAX=0.0 
DO 5 1=1,961 
K=I-1 

IF (1-1) 2,2,4 

2 GG=1.0 
GO TO 1 

4 GG=ABS (SIN (PI*K*416, 67/400000) /(PI*K*416. 67/400000) ) 

G=GG 

1 CONTINUE 

C=COS (2*PI*K*DEL) 

S=SIN (2*PI*K«DEL) 

CC=COS (4*PI*K*DEL) 

SS=SIN (4*PI*K*DEL) 

XR=A0+A1*C+A2*CC 

XI=A1*S+A2*SS 

Xa=1+B1*C+B2*CC 

YI=B1*S+B2*SS 

R=XR*YE+XI*YI 

RI=XB*YI-XI*YR 

D= YR**2+YI**2 

H (I)= (SQRT (R**2+RI<'*2) /D) *GG 
IF (H(I)-HMAX) 6,3,0 

3 HMAX=H(I) 
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6 CONTINDE 

5 CONTINOE 

DO 9 1=1,961 
K=I-1 

F (I) =K=*0.4 
A (I) =H (I) /HMAX 
9 CONTINDE 

CALL INIT 
CALL XFEM(4) 
CALL YFEH(4) 
CALL NETS (961) 
CALL PLOT(?,A) 
CALL PAUSE 
CALL FIN 
END 
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APPENDIX C 



EXPERIMENrAL DATA FOR SNR GRAPHS 



Noise Volts 
(RMS) 



Signal Volts 
(RMS) 



Coef f . 
b 

1 

+ 0.8 
+ 0.7 
+ 0.6 
+ 0.5 
+ 0.2 
0 . 0 
- 0.2 
-0.5 
- 0.6 
-0.7 
- 0.8 



V 

o 

0.061 
0.059 
0.056 
0.056 
0.055 
0.063 
0.068 
0.079 
0.097 
0.118 
0. 1 14 



V 

i 

0.038 
0.405 
0.375 
0.039 
0.037 
0.039 
0. 037 
0. 039 
0.038 
0.042 
0. 03 8 



V 

o 

0. 036 
0.035 
0. 035 
0.036 
0.039 
0. 039 
0.043 
0.062 
0. 083 
0. 114 
0.195 



V 

i 

0.021 

0.024 

0.202 

0.226 

0.023 

0.022 

0.020 

0.022 

0.021 

0.023 

0.021 



In 

-5.36 
-4.51 
-5.37 
-4.76 
-3.90 
-5. 18 
-5.37 
-5.07 
-5.25 
-5.21 
-5.36 



SNR 

(DB) 

Out 

-4.58 
-4.55 
-4.00 
-3.84 
-2. 99 
-4. 17 
-3.98 
-2. 1 3 
-1.31 
-0.30 
+ 2.63 



*Data plotted as dots with squares around them 
in Figure 5.26. 



Improv. 

+0.78 
-0.04 
+ 1.37 
+0.92 
+0.91 
+ 1 .02 
+ 1.39 
+2.95 
+ 3.94 
+4.91 
+7.99 
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Noise Volts 
(RMS) 



Signal Volts 
(RMS) 



Coeff . 
b 

1 

+ 0.8 
+0.7 
+ 0.6 
+ 0.5 
+ 0.4 
+0.3 
+ 0.2 
+ 0.1 
0.0 
- 0.1 
- 0.2 
-0.3 
-0.4 
-0.5 
- 0.6 
-0.7 
- 0.8 
-0.9 



V 

o 

0.126 
0.115 
0.119 
0.120 
0.123 
0.125 
0.125 
0.135 
0.134 
0. 135 
0.150 
0.152 
0.162 
0. 196 
0.219 
0.244 
0.320 
0.575 



V 

i 

0.082 
0. 082 
0. 085 
0.085 
0. 085 
0. 085 
0.085 
0. 085 
0. 085 
0. 085 
0. 085 
0. 085 
0.085 
0. 085 
0. 085 
0. 082 
0.082 
0. 085 



0. 066 
0. 060 
0.060 
0.063 
0.070 
0.077 
0.076 
0.088 
0.094 
0. 088 
0.099 
0. 098 
0. 125 
0.181 
0.218 
0.295 
0.425 
1.300 



0.044 
0.044 
0.044 
0.045 
0.044 
0.045 
0.045 
09045 
0.045 
0.045 
0.045 
0.045 
0.044 
0.045 
0.044 
0. 044 
0.044 
0.044 



In 

-5.53 

-5.43 

-5.72 

-5.58 

-5.70 

-5.60 

-5.62 

-5.60 

-5.62 

-5.62 

-5.62 

-5.60 

-5.70 

-5.58 

-5.72 

-5.43 

-5.53 

-5.72 



SNR 

(DB) 

Oat 

-5.68 
-5.67 
-5. 96 
-5.61 
-4.91 
-4.26 
-4.38 
-3.77 
-3.08 
-3.72 
-3.61 
-3.81 
-2. 25 
-0.69 
-0.04 
+ 1. 65 
+ 2.46 
+ 7.09 



♦Data plotted as dots with triangles around thea 
in Figure 5.26. 



laprov. 

-0.15 
-0.24 
-0.24 
-0.0 3 
+0.79 
+ 1.34 
+ 1 .24 
+ 1.84 
+2.54 
+ 1.90 
+ 2.01 
+ 1.79 
+3.45 
+4 . 89 
+5.68 
+7.08 
+7.99 
+ 12. 23 
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Noise Volts 


Signal Volts 




SNR 


(RMS) 


(RMS) 




(DB) 


V V 

o i 


V V 

o i 


In 


Out Iraprov 



Coeff . 
b 

1 

+ 0.7 
+ 0.6 
+ 0.5 
+ 0.4 
+ 0.3 
+ 0.2 
+ 0.1 
0.0 
- 0.1 
- 0.2 
-0.3 
-0.4 
-0.5 
- 0.6 
-0.7 
- 0.8 
-0.9 



0.204 
0.206 
0.207 
0.212 
0.220 
0.228 
0.236 
0.235 
0.236 
0.24 2 
0.26 1 
0.280 
0.312 
0.353 
0.41 0 
0.525 
0.905 



0. 150 
0. 150 
0.150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 
0. 150 



0. 038 
0. 037 
0.037 
0. 039 
0. 042 
0.043 
0. 047 
0.045 
0.047 
0.049 
0. 056 
0. 063 
0. 073 
0. 094 
0.123 
0. 1 97 
0.480 



0.023 
0.023 
0.023 
0.023 
0.023 
0. 023 
0.023 
0.023 
0.023 
0.023 
0.023 
0.023 
0.023 
0.023 
0.023 
0.023 
0.023 



-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 

-16.5 



-14.6 

-14.9 

-14.9 

-14.6 

-14.4 

- 14.5 

-14.0 

-14.3 

-14.0 

-13.9 

-13.4 

-12.9 

-12.7 

-11.5 

-10.5 

-8.51 

-5.51 



+ 1.86 
+ 1.56 
+ 1 . 59 
+ 1.93 
+2.09 
+ 1.99 
+2.46 
+ 2. 18 
+2.46 
+ 2.61 
+ 3.1 1 
+ 3.52 
+3.80 
+4.98 
+ 6 .02 
+7.96 
+ 10. 9 



♦Data plotted as dots with circles around them 
in Figure 5.26. 
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